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ABSTRACT 


Gold  was  found  to  diCfuso  into  silicon  by  a  complox  mocbanism 
involving  a  vacancy  •controlled  interstitial-substitutional  equilibrium. 
This  led  to  very  complex  diffusion  concentration  profiles.  In  analysing 
the  problem,  a  new  experimental  value  was  found  for  the  self-diffusion 
coefficient  of  silicon,  1.81  X  10^  exp  (-112^:20  kcal/RT)cm^/sec,  which 
compared  very  favorably  with  previous  data  on  diffusion  of  bismuth, 
germanium,  and  tin  in  silicon.  The  interstitial  gold -diffusion  coeffi¬ 
cient  was  found  to  be  2.4X10*^  exp  (-8.9^2  kcal/RT),  the  substitutional 
gold -diffusion  coefficient  to  be  2.  75  X  10"^  exp  (-47  ±  10  kcal/RT),  the 
equilibrium-interstitial  gold  solubility  to  be  5.95  X  10^^  exp  (-58  O  10 
kcal/RT),  and  the  equilibrium-substitutional  gold  solubility  (Collin’s 
d«te)  below  1200°C  to  be  8.  15  X  10^^  exp  (-40.6  kc»l/RT). 
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I.  INTRODUCTION 


Gold  conveniently  furnishes  recombination  centers  in  silicon. 

Because  of  this,  gold  has  been  deliberately  diffused  into  silicon  to 

make  high-speed  semiconductor  devices.  The  diffusion  of  gold  into 

silicon  was  first  studied  in  detail  by  Struthers  using  radiotracers.  ^ 

The  diffusion  coefficients  measured  gave  a  good  Arrhenius  plot  from 

800^C  to  IZOO^C  to  yield  D  =  0.011  exp  (-25,  800/RT).  More  recent 

work,  however,  has  indicated  that  gold  diffuses  into  silicon  at  a  much 
7  8 

slower  rate.  '  Limited  evidence  that  the  diffusion  may  be  structure- 

7 

sensitive  was  also  reported.  It  was  observed  by  Dash  that  gold  diffu¬ 
sion  at  1150^0  resulted  in  the  climb  of  dislocations.  ^  Dash  suggested 
that  this  phenomenon  could  be  explained  by  the  substitutional- interstitial 
diffusion  mechanism  hypothesized  for  diffusion  of  copper  into 
germanium.  In  this  mechanism,  the  diffusing  atoms  move  in  two 

manners:  interstitially  and  substitutionally.  A  high  interstitial  diffu¬ 
sion  rate  is  combined  with  low  interstitial  solubility.  Conversely,  a 
low  substitutional  diffusion  rate  is  combined  with  high  substitutional 
solubility.  The  net  result  is  that  the  diffusive  flux  is  carried  almost 
entirely  by  interstitial  atoms.  The  interstitials  then  combine  with 
vacancies  to  form  substitutional  atoms.  The  consiiming  of  vacancies 
leads  to  dislocation  climb. 
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The  purpose  of  the  work  described  in  this  paper  was  to  determine, 
quantitatively,  the  mechanism  of  gold  diffusion  into  silicon.  Measure¬ 
ments  of  the  diffusion  were  carried  out  from  700^C  to  1300^C  using 
198  199 

radiotracer  Au  and  also  Au  .  Almost  100  wafers  were  diffused, 

including  dislocation- free  Czochralski  silicon  and  floating- zone 

4  2 

material  (having  about  10  dislocations/cm  ),  high  and  low  resistivity 

13 

n-  and  p-type  silicon,  and  previously- diffusion- doped  samples. 

The  substitutional- interstitial  mechanism  was  verified,  and  the  param¬ 
eters  involved  were  determined.  Values  were  determined  for  inter¬ 
stitial  gold  solubility,  interstitial  gold  diffusion  coefficient,  substitu¬ 
tional  gold  diffusion  coefficient,  and  the  self- diffusion  rate  of  silicon. 

II.  EXPERIMENTAL  PROCEDURE 

The  experimental  techniques  used  in  this  work  are  reported  else- 

14 

where  in  much  greater  detail  than  is  attempted  here.  Both  (111) 

4  2 

floating- zone  grown  silicon  of  about  10  dislocations /cm  and  dislocation- 
free  (100)  Czochralski- grown  silicon  were  used.  Comparison  with 
dislocation- free  (111)  Czochralski- grown  material  showed  that 
crystal  orientation  did  not  affect  the  gold  diffusion.  Therefore,  in  the 
figures  to  follow,  (100)  is  taken  to  mean  Czochralski  material  and  (111) 
to  mean  floating- zone  material.  The  specimens  were  in  the  form  of 
thick  cylindrical  wafers  of  about  1  to  2  cm  in  diameter  and  0.  5  cm  in 
thickness.  The  wafers  were  parallel-lapped  with  13-micron  garnet.  Extreme 
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care  was  taken  to  attain  a  very  clean,  oxide -free  surface  before  plating 

198  199 

of  the  gold  which  was  tagged  with  radiotracers.  Both  Au  and  Au 

were  used  at  different  times  in  tagging  stable  gold -chloride  solutions 

199 

with  radiotracer.  Nearly  carrier -free  Au  was  also  used,  undiluted 
(as  supplied  by  Oak  Ridge  National  Laboratory),  in  some  experiments 
to  produce  very  low  surface  concentrations  of  gold.  The  gold  was  chemi- 
plated  onto  the  wafers  by  injecting  a  dilute  gold  chloride  solution  into 
hydrofluoric  acid  on  the  surface  of  the  wafer.  Autoradiographs  before 
and  after  diffusion  showed  that  a  very  uniform  layer  was  obtained  (Ref. 

14).  The  wafers  were  then  encapsulated  in  quartz  in  an  argon  atmosphere 
and  placed  in  the  diffusion  furnace  controlled  to  ±  2^C  for  a  period  of  time 
ranging  from  15  minutes  to  several  days.  After  diffusion,  the  capsules 
were  quenched  by  placing  them  on  a  massive  brass  block.  Excess  sur¬ 
face  activity  was  removed  by  dissolving  the  gold  with  aqua  regia.  Count¬ 
ing  was  then  begun  using  a  Geiger -Muller  tube.  Counting  accuracy  was 
usually  on  the  order  of  a  few  percent.  A  collimator  on  the  wafer  was  used 
to  eliminate  edge  effects  and  non-parallelism.  Thin  layers  were  removed 
using  a  high  precision  machine  with  fine  silicon-carbide  paper.  Thick¬ 
nesses  were  measured  by  a  Mikrokator  and  by  a  Mettler  microbalance  to 
-5 

about  5X10  cm.  Best  results  were  obtained  by  counting  the  wafers 
rather  than  the  paper.  The  counting  rate  versus  distance  was  converted 
to  a  concentration  profile,  using  the  techniques  described  in  Ref,  14. 
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HI.  EXPERIMENTAL  RESULTS 


The  net  results  of  the  experiments  were  plots  of  gold  concentration 
versus  distance  into  the  silicon.  The  initial  experiments  were  performed 
between  SOO^C  and  IIOO^C  on  dislocation-free  silicon.  The  diffusion  was 
much  slower  than  reported  by  Struthers,^  as  shown  in  Fig.  1.  In  addition, 
the  concentration  profiles  did  not  remotely  resemble  a  complementary- 
error  function  distribution,  which  would  be  expected  for  a  saturated  sur¬ 
face  condition  with  a  constant  diffusion  coefficient.  It  was  found,  however, 
that  when  plotted  on  log-log  paper  the  profiles  were  all  nearly  identical. 

An  empirical  curve  of  the  results  is  shown  in  Fig.  2.  An  interesting  ob¬ 
servation  is  that  the  portion  of  this  curve  from  a  concentration  of  about 
10  to  98  percent  of  the  solubility  gives  a  very  good  straight-line  plot  of 
log  concentration  versusv^c/Kt.  The  values  of  K  for  various  materials  and 
diffusion  conditions  are  given  in  greater  detail  in  Ref.  14.  The  effect  of 
high  prior  doping  with  n-  and  p-type  impurities  is  also  discussed  in  Ref.  14. 
Phosphorous -diffused  layers  with  high  surface  concentrations  were  found 
to  both  hinder  the  gold  diffusion  rates  and  to  increase  the  gold  solubilities. 

Further  experiments  revealed  that  the  diffusion  rates  were  higher 
when  floating -zone  silicon  was  used  than  when  dislocation-free  silicon  was 
used.  It  was  also  found  that  the  ideal  error -function  complement  profile 
was  obtained  for  high  temperatures  (1300°C),  Fig.  3,  and  for  very  long  dif¬ 
fusion  times  at  lower  temperatures  (700°C),  Fig.  4.  In  Fig,  3,  note  the 

-3 

evidence  of  out -diffusion  within  10  cm  of  the  specimen  surface.  This 
out -diffusion  is  caused  by  reduced  solubility  during  cooling  from  the 
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diffusion  temperature  at  the  end  of  the  run.  As  would  be  expected,  the  out- 
diffusion  was  more  severe  in  the  high-dislocation-density  silicon.  Figure  3 
also  shows  Gaussian  distributions  obtained  under  conditions  of  very  low 
surface  concentration,  designated  C.  F.  .n  the  figures  (for  nearly  carrier- 
free  radiotracer).  It  was  possible  to  obtain  Gaussian  distributions  at 
1200^C>  Fig,  5,  and  even  at  lOOO^C  under  conditions  of  very  low  surface 
concentrations  and  high  dislocation  density.  Figure  5  also  shows  that  the 
diffusion  profiles  for  dislocation-free  silicon  had  already  deviated  from 
ideal  and  were  approaching  those  of  Fig.  2.  Figure  6  shows  typical  con¬ 
centration  profiles  obtained  for  diffusion  at  llOO^C.  Note  that  even  the 
profiles  for  the  (111)  floating-zone  silicon  deviated  from  the  ideal  error- 
function  complement  at  llOO^C. 

From  the  foregoing  results,  diffusion  coefficients  were  calculated  and 
are  shown  in  Fig.  7  as  log  D  versus  1/T  (an  Arrhenius  plot).  Calculation 
of  diffusion  coefficients  for  the  error-function  complement  (erfc)  and 
Gaussian -shaped  profiles  were  easily  obtained  by  superimposing  log -log 
plots  of  these  curves  on  the  data.  The  more  complex  experimental  profiles 
were  analyzed  by  similarly  superimposing  on  a  log-log  erfc  plot  those  data 
points  which  seemed  to  fit.  Many  of  the  profiles,  for  example,  had  erfc 
tails,  as  illustrated  in  Fig.  6.  The  erfc  tails  at  800^C  and  900°C  were 
much  more  pronounced  than  in  Fig.  6,  Extrapolation  of  the  erfc  fit  to  the 
tails  back  to  the  surface  of  the  sample  also  gave  values  for  a  solubility. 

As  will  be  seen  later,  these  tails  are  believed  to  be  caused  by  interstitial 
diffusion  alone.  Hence,  the  extrapolated  surface  concentration  is  believed 
to  be  the  interstitial  gold  solubility.  The  ratio  of  total  gold  solubility  to 
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the  interstitial  gold  solubility  is  given  in  Fig.  8j  where  it  is  compared 
with  values  obtained  from  calculations  using  data  in  Fig.  7.  Because 
the  interstitial  solubility  is  so  much  lower  than  the  substitutional  solu¬ 
bility,  the  substitutional  solubility  is  very  nearly  equal  to  the  total  gold 
solubility.  Fitting  the  erfc  curve  to  the  rest  of  the  data  of  complex  pro¬ 
files  likewise  gave  diffusion  coefficients,  shown  in  Fig.  7.  The  diffusion 
coefficients  given  by  Boltaks  '  and  Struthers  are  also  shown  in  Fig.  7. 

IV.  DISCUSSION 

In  attempting  to  determine  the  best  mechanism  for  gold  diffusion  in 
silicon  based  on  the  experimental  data,  we  will  assume  a  substitutional- 
interstitial  process  and  analyze  the  data  on  this  basis.  If  the  results 
are  reasonable,  this  can  be  taken  as  strong  evidence  supporting  the 
assumed  mechanism. 

The  various  types  of  diffusion  coefficients  which  can  arise  with  the 
assumed  mechanism  are  summarized  in  Table  1.  The  method  of  deri¬ 
vation  is  illustrated  in  the  Appendix. 
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Table  1.  Diffusion  Coefficients  Arising  in  Substitutional -Interstitial 

Diffusion  Mechanisnr) 


Type  of  Curve 

Case 

Assumptions^ 

Diffusion  Coefficient^ 

Erfc  Overall 

(i) 

C  =  {C  )  .  high  N,. 
v  '  v's  ®  dis 

and/or  high  T 

(2) 

Long  time,  low  T 

°sub 

(3) 

C.  =  {C.)  ,  low  N,. 
i  'is  dis 

(C  ) 

V'S  rx 

|C  )  4-  (C  ■J'  V 

'  v's  sub  « 

Gaussian 

(4) 

Low  surface 

concentration,  C.  «  C 

1  V 

ic.)  +IC'  J 

'  I's  sub  s 

(5) 

low  surface 

concentration,  C.  »  C 

1  V 

D. 

1 

Tail  Region 

(6) 

C.  «  C 

1  V 

<^i>. 

(C.)  +(C  .) 

'  I's  sub's 

(7) 

C.  »  C 

1  V 

D. 

1 

Surface  Region 

(8) 

Long  time,  low  T 

^sub 

(9) 

(Cy), 

(C  )  +TC  J 
'v's  '  sub's 

denotes  dislocation  density,  T  denotes  temperature,  C  denotes 
concentration,  D  denotes  diffusion  coefficient.  Subscripts  as  follows; 
i  =  interstitial  gold,  sub  =  substitutional  gold,  v  =  vacancy,  s  =  satura¬ 
tion  or  solubility  value. 
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If  the  diffusion  coefficients  and  concentrations  of  vacancies, 
interstitial  gold,  and  substitutional  gold  can  all  be  assumed  to  vary 
according  to  the  Arrhenius  law,  then  a  plot  of  log  D  versus  1/T 
should  show  a  straight  line.  Figure  7  shows,  in  fact,  three  apparent 
straight  lines,  which  are  designated  as  Curves  I,  II,  and  III.  The  equa* 
tions  of  these  lines  in  cm^/sec  are 


Curve  I:  D  =  2.  75  X  10"^  exp  (-47,  400/RT)  ,  (1) 

Curve  II:  D  =  1.  15  X  10^  exp  (-71, 800/RT)  ,  (2) 

Curve  III:  D  =  1.  78  X  10'^  exp  (-26,  000/RT)  .  (3) 


The  problem  is  to  determine  the  meaning  of  these  curves.  From 
Table  1,  it  is  observed  that  four  different  diffusion  coefficients  might 
arise.  We  must  determine  which  of  the  coefficients  is  applicable  to 
each  measured  curve. 

If  we  consider  the  vacancy^  diffusion- controlled  situation  (Case  3, 
Table  1),  then  the  effective  diffusion  coefficient  is  given  by 


(C  ) 

'  sub' s 


(4) 
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Here,  it  is  assumed  that  (C^)^  «  which  is  reasonable.  Note, 

however,  that 


°v  =  'V.  (^)  D. 

■  <5) 

"  °sd  ('W) 

where  (X^)^  is  the  equilibrium  atom  fraction  of  vacancies  in  silicon. 

N  is  Avogadro's  number,  p  is  the  density  of  silicon,  and  W  is  the 
atomic  weight  of  silicon.  The  term  (Np/W)  is  actually  the  number  of 
silicon  atoms  per  unit  volume. 

Substituting  Eq.  (5)  into  Eq.  (4)  and  solving  for  D  ,  we  obtain 

sd 

°sd  "  sub's  °Au  •  (6) 

If  it  is  assumed  that  in  fact  shown  in  Fig.  8, 

then  the  equilibrium  substitutional  gold  solubility  is  very  nearly 

equal  to  the  total  gold  solubility,  as  measured  by  Collins  et  al.  ^  and  by 
Struthers.  With  this,  I^g^  can  be  calculated  using  Eq.  (6)  and  the 
curves  I,  II,  and  III,  in  turn.  The  true  result  should  be  approximately 
equal  to  the  diffusion  coefficients  known  for  substitutional  elements  in 
silicon,  and,  in  particular,  for  the  Group  IV  elements,  germanium  and 
tin. 
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Collins'  data  for  gold  solubility  can  be  given  by 

(C8ub)g  =  8.  15X  10^^  exp  (-40,  600/RT)  (gold  atoms/cm^  Si)  .  (7) 

This  is  particularly  accurate  below  IZOO^C. 

Substituting  Eq.  (7)  and  Eq.  (1),  (2),  and  (3),  in  turn,  into  Eq.  (6), 
we  obtain,  respectively: 


Curve  I: 

D  ,  =  4.  32  X  10"^  exp  (-88,  OOO/RT)  , 

SQ 

(8) 

Curve  II: 

=  1.81  X  10^  exp  (-112,400/RT)  , 

(9) 

Curve  III; 

=  2.  8  X  10'^  exp  (-67,  200/RT) 

(10) 

As  shown  in  Fig.  9i  the  silicon  self- diffusion  coefficient  deter¬ 
mined  from  Curve  lli  Eq.  (9)i  is  very  close  to  the  diffusion  coefficients 
of  bismuth,  tin,  germanium,  and  fairly  near  that  for 

indium^  ^  in  silicon. 

The  gold  solubility  given  by  Struthers^  differs  from  the  data  of 
Collins*  used  above.  Struthers  found,  approximately. 
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This  gives,  for  Curve  II: 


Dsd  =  6.  63  X  10  exp  (-87,  100/RT)  .  (12) 

This  is  also  shown  in  Fig.  9<  It  appears  not  to  give  as  good  an  agree¬ 
ment  as  that  found  using  the  solubility  data  of  Collins.  Although  this 
is  not  entirely  certain,  we  shall  tentatively  conclude  from  this  that 
Collins'  data  for  gold  solubility  is  better  than  that  of  Struthers. 

On  the  basis  of  the  above  calculations  and  comparisons,  we  con¬ 
clude  that  Curve  II  of  Fig.  7  represents  the  vacancy- controlled  diffu¬ 
sion 


D 


The  only  other  possible  explanation  is  that  this  curve  is,  instead, 


D. 

1 


In  such  an  event,  Curve  III  would  have  to  represent  D..  The  ratio  of 
values  on  Curve  III  to  those  on  Curve  II  must  then  yield  values  for 
^^sub^s^^^i^s  true.  These  values  are  to  be  compared  with 

the  experimental  ratios  given  in  Fig.  8.  Because  there  is  virtually 
no  agreement, we  have  shown  that  Curve  II  cannot  represent  the 
interstitial-controlled  diffusion.  Therefore  Curve  II  must  represent 


13 


D 


V 


as  originally  concluded.  Curve  III  must  be 


(C.) 

I's 

(cr""7r 

'  sub' s 


D. 

1 


It  cannot  be  D.  because,  as  we  have  just  shown,  this  would  necessitate 
a  lower  curve  for 


which  is  not  present.  It  cannot  be  because  the  values  are  far  too 

large  for  substitutional  diffusion  in  silicon.  Using  the  experimental 
results  in  Fig.  8  for  (C  ,  )  /(C.)  (as  determined  from  the  tails  on 
the  concentration  profiles)  together  with  Eq.  (3)  for  Curve  III,  we  find 
that 


D.  =  2.  44  X  10"'^  exp  (-8,  900/RT)  (13) 

From  this  equation  we  calculate  that,  at  1200°C,  D.  =  1. 3X  10"^  and, 
at  1000°C,  D.  =  7.  5  X  lO"^  . 
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This  compares  very  well  with  the  interstitial  diffusants:  lithium, 
nickel,  copper,  helium,  and  hydrogen^^  in  silicon. 

Now  only  Curve  I  remains.  By  the  foregoing  process  of  elimina¬ 
tion  and  by  virtue  of  the  low  values,  Curve  I  must  represent  the 

pure  substitutional  diffusion. 

As  a  final  point,  it  is  worthwhile  to  re-examine  the  older  data  for 

gold  diffusion  in  silicon.  Struthers*  data^  are  shown  in  Fig.  7  as  the 

7 

S- points  and  Boltaks'  by  B- points.  Struthers*  data  fall  very  nearly 
on  Curve  III  for  interstitial- controlled  diffusion.  From  this,  it  is  sur¬ 
mised  that  the  original  experimental  conditions  were  probably  as 

3 

follows:  very  low  perfection  silicon  with  more  than  10  dislocations/ 

2 

cm  and/or  a  surface  concentration  below  saturation.  Struthers  has, 

Z  5 

in  fact,  verified  that  both  of  these  were  his  experimental  conditions. 
Boltaks*  data,  on  the  other  hand,  fall  very  nearly  on  the  vacancy- 
controlled- diffusion  Curve  II,  indicating  that  he  used  saturated  surface 
conditions. 


V.  CONCLUSIONS 

It  has  been  determined  that  gold  diffuses  into  silicon  by  a  mecha¬ 
nism  analogous  to  that  of  copper  diffusing  into  germanium.  Movement 
of  gold  atoms  primarily  involves  an  interstitial  diffusion,  controlled 
by  the  concentration  of  vacancies  (which  are  necessary  to  attain  the 
higher  substitutional  concentrations)*  Such  a  complex  mechanism  often 
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leads  to  unusual  diffusion- concentration  profiles*  as  exemplified  by 
the  most  common  one,  shown  in  Fig.  2.  Prediction  of  the  exact  con¬ 
centration  profile,  a  priori,  would  necessitate  a  numerical- computer 
solution  of  the  differential  equations  describing  simultaneous  diffusion 
of  interstitials,  substitutionals,  and  vacancies.  Since  this  mechanism 
usually  causes  a  strong  dependence  of  the  diffusion  rate  on  crystalline 
perfection  and  on  quantity  of  diffusant  the  aforementioned  computer 
solution  for  a  diffusion- concentration  profile  would  necessitate  knowl¬ 
edge  of  the  rate  of  vacancy  generation  by  dislocations  in  silicon. 

Values  for  the  latter  have  not  yet  been  determined. 

In  the  course  of  these  analyses,  the  parameters  of  fundamental 
importance  in  the  diffusion  of  gold  into  silicon  have  been  determined 
as  follows. 

2 

Self- diffusion  coefficient  of  silicon  (cm  /sec).  Fig.  4: 

=  1.  81  X  10^  exp  (-112  ±  20  kcal/RT)  (9) 

2 

Interstitial  gold- diffusion  rate  in  silicon  (cm  /sec): 

D.  =  2.  44  X  10"^  exp  (-8.9  ±  2  kcal/RT)  (13) 

2 

Substitutional  gold- diffusion  rate  in  silicon  (cm  /sec): 

=  2.  75X  10"^  exp  (-47  ±  10  kcal/RT)  (1) 
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I  3 

Eqxiilibrium- subatitutional  gold  solubility  ( atoms /cm  )  below 
1200°C: 

(C^^b^s  =  ®-  ^  kcal/RT)  (7) 

3 

Equilibrium- interstitial  gold  solubility  (atoms /cm  ): 

(C.)^  =  5.  95  X  10^“^  exp  (-58  ±  10  kcal/RT)  (14) 
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APPENDIX:  TYPICAL  DIFFUSION  SOLUTIONS 


Following  are  two  typical  solutions  for  limiting  cases  of  a  vacancy- 
controlled  interstitial- substitutional  diffusion.  The  reaction  of  an  in¬ 
terstitial  gold  atom  going  into  a  substitutional  position  is  represented 
by 


Au.  +  V 


Au 


sub 


(15) 


where  V  is  a  vacancy. 

It  is  assumed  that  the  concentration  of  these  three  is  at  equilibrium 
at  all  times.  This  is  expressed  by 


^sub 

TTjr 

X  V 


=  K 


(16) 


where  K  is  an  equilibrium  constant  which  is  a  function  of  temperature 
alone. 


Case  A 


Assume  C^  =  (C^)^,  the  maximum  solubility  of  vacancies: 


D  ,  «  D  ;  and  D  ,  «  D. 
sub  V  sub  1 
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4 


The  first  condition  holds  for  a  very  high  dislocation  density  (since 
vacancies  are  produced  at  dislocations)  or  where  the  gold  concentra¬ 
tion  is  so  low  that  few  vacancies  are  consumed.  The  gold  concentra¬ 
tion  is  also  very  low  for  the  first  gold  atoms  that  penetrate  deeply  into 
the  crystal.  The  first  condition  may  also  be  true  even  in  dislocation- 
free  crystals  at  high  temperatures  where  the  number  and  rate  of  gener¬ 
ation  of  vacancies  is  greater.  The  second  and  third  conditions,  which 
state  that  the  substitutional  atoms  are  relatively  stationary,  are  known 
to  be  at  least  approximately  true. 

The  expression  for  interstitial  diffusion  is 


8C. 

1 

at 


8C 


sub 

'St— 


(17) 


But,  from  Eq.  (l6),  we  know  that 

C  ,  =  KC.  (C  )  (18) 

sub  1  '  v's 

Therefore,  Eq.  (17)  becomes 


8^C.  r  -|  8C. 


(19) 


For  constant  interstitial  surface  concentration,  the  solution  to  this  is 


C.  =  (C.)  erfc 
1  '  I's 


2[(D.)^,,t] 


m 


(20) 
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where 


D. 


(D.) 


'i'eff  “  I  +K(C) 

V  s 


(C.) 

'  I's 


(C.)  +(C  ”5  °i  • 

'  I's  '  sub's 


(21) 


Using  the  equilibrium  expression,  Eq.  (18),  we  find  that 


C  u  =  (C  ,  )  erfc 
sub  '  sub' s 


m 


(22) 


Hence,  the  effective  diffusion  coefficient  for  gold  is 


(C;). 


(C;). 


D:  . 


(23) 


assuming 


(C.)  «  (C  J 
'is  sub's 


Case  B 


Assume  D.  »  D  ,  ;  then 
1  sub 


D. 

1 


ac.  8C 


and  that  the  silicon  is  dislocation  free. 


(24) 
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In  essence,  it  is  assumed  that  only  the  interstitial  gold  atoms 
move,  so  that  the  vacancy  concentration  is  the  saturated  value  reduced 
by  the  number  of  substitutional  gold  atoms,  i.  e.  , 


C 


v 


(C  ) 

'  v's 


sub 


(25) 


Substituting  this  into  Eq.  (6)  we  find  that 


'sub 


=  K 


(26) 


Differentiating  with  respect  to  time,  there  results 

»=.»b  C.(C^).-]  8C; 

-5t—  [1  +KC,  -  ,i  tKC.H 


(27) 


Substituting  this  into  Eq.  (17),  we  obtain 


8^C.. 


ax‘ 


1  8Ci 


(28) 


where 


(1  +KC.)^ 

f(C.)  =  - - 

^  (1  +KC.)'' +  K  (C^)^ 


(29) 
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I 


For  constant  interstitial  surface  concentration,  the  solution  to  this  is 


i 


c. 

1 


(C.) 

IS 


gj - — 


m 


26 

(30) 


where 


ThiSi  in  turn,  gives  an  effective  diffusion  coefficient  for  the  total  gold 
content  which  is  a  function  only  of  concentration.  This  is  not  much 
help  quantitatively  but  indicates  that  C/C^  can  be  a  single  function  of 
x/\[i  at  a  given  temperature  (as  was  found  experimentally  in  Fig.  2) 
and  still  conform  to  the  present  mechanism. 
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Fig.  1.  Concentration  Profile  for  Gold  Diffusion 
into  Dislocation-Free  Silicon  at  1097®C; 
Compared  to  Struthers'  Results 
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mUTIVi  MiTANOi,  i/yRT 


Error -function  complement  curve  is  shown  for  comparison 
of  shape  only,  relative  position  being  unimportant.  Here 
C  represents  the  gold  concentration  at  distance  x  from  the 
surface  (where  the  concentration  is  equal  to  the  solubility 
value  Cg),.  t  is  the  time  of  diffusion  and  K  is  a  constant 
necessary  to  superimpose  the  data,  and  maybe  considered 
to  be  an  effective  diffusion  coefficient. 


Fig.  2.  Empirical  Correlation  of  Concentration  for  Gold 
Diffusion  into  Dislocation-Free  Silicon 
from  800OC  to  llOO^C 
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199 

C.  F.  denotes  that  undiluted  radiotracer  (Au  as  supplied)  was  used 
in  plating  surface  before  diffusion. 


Fig.  3.  Concentration  Profiles  for  Gold  Diffused  into  Silicon  at  1300°C 
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Fig.  5.  Distribution  of  Gold  Diffused  into  Silicon  at  1200°C 
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Distribution  of  Gold  Diffused  into  Silicon  at  1000  C 
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Fig,  7.  Diffusion  Coefficients  Calculated  from 
Experimental  Diffusion  Concentration 
Profiles  as  a  Function  of  Inverse 
Temperature  (Arrhenius  Plot) 
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INTERSTITIAL  GOLD  SOLUBILITYCC,^ /(C;), 
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Fig.  8.  Ratio  of  Substitutional  to  Interstitial  Gold  Concentrations,  as 
Calculated  from  Experimental  Diffusion  Concentration 
Profiles  and  from  Curves  II  and  III  in  Fig.  7 
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Fig.  9.  Comparison  of  Calculatod  Salf-Diffuaion 
Coafflciant  of  Silicon  with  Pravlous  Data 
for  Impurity  Diffusion  into  Silicon 
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